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Laser interferometric gravitational-wave interferometers implement Fabry-Perot cavities to increase their 
peak sensitivity. However, this is at cost of reducing their detection bandwidth, which origins from the propa¬ 
gation phase delay of the light. The “white-light-cavity” idea, first proposed by Wicht et al. [Optics Commu¬ 
nications 134 , 431 (1997)], is to circumvent this limitation by introducing anomalous dispersion, using double- 
pumped gain medium, to compensate for such phase delay. In this article, starting from the Hamiltonian of 
atom-light interaction, we apply the input-output formalism to evaluate the quantum noise of the system. We 
find that apart from the additional noise associated with the parametric amplification process noticed by others, 
the stability condition for the entire system poses an additional constraint. Through surveying the parameter 
regimes where the gain medium remains stable (not lasing) and stationary, we find that there is no net enhance¬ 
ment of the shot-noise limited sensitivity. Therefore, other gain mediums or different parameter regimes shall 
be explored for realizing the white light cavity. 

PACS numbers: 


I. INTRODUCTION 


Second-generation large-scale interferometric gravitational 
wave (GW) detectors, such as advanced LIGO m, advanced 
VIRGO d and KAGRA (3], are designed to operate at better 
sensitivity than the first generation detectors. This improve¬ 
ment in sensitivity comes from increase in the optical power 
and introduction of a signal recycling mirror (SRM) to the 
initial configuration @. The SRM on the dark port forms a 
signal recycling cavity (SRC) with the input test mass mir¬ 
ror (ITM). The position of SRM determines the propagation 
phase of the signal light inside the SRC, and control of the 
SRM parameters allows for adjustments to the frequency re¬ 
sponse of the interferometer mm. Two typical operational 
modes are the signal recycling mode and resonant sideband 
extraction mode (RSE). The signal recycling mode enhances 
the sideband carrier of the GW signal inside the cavity, while 
the RSE mode increases the detection bandwidth which is the 
effective bandwidth of the combined SRC and arm cavity (7). 

However, broadening the detection bandwidth in the RSE 
mode comes at the loss of the peak sensitivity; while enhanc¬ 
ing the peak sensitivity in SRC results in a narrower detection 
bandwidth. This trad-off is represented by the integrated sen¬ 
sitivity: 
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which only depends on the intra-cavity power P c and cavity 
length L arm , and is independent of the property of SRC. The 
(Do, 0 )fsr 5 c here are the laser frequency, free-spectral range, 
and the speed of light, respectively. Here, we only consider 
the shot noise limited strain sensitivity 5/^ (£2) since radia¬ 
tion pressure noise can in principle be evaded using frequency 
dependent readout or sufficiently heavy test masses. Such a 
trade-off between bandwidth and peak sensitivity is due to 
the accumulated phase of the sideband field propagating in¬ 
side the arm cavity. There are several proposals in the litera¬ 
ture that try to achieve the resonant amplification of the sig¬ 


nal without decreasing the bandwidth, using the idea of white 
light cavity. Among those, Wicht et al. were the first to sug¬ 
gest placing an atomic gain medium with anomalous disper¬ 
sion inside the SRC to cancel the propagation phase El. 
This idea was then followed by Pati and Yum et al. with dif¬ 
ferent types of active mediums iflOlITTll . 

The anomalous dispersion phenomenon and the interesting 
“superluminal” physics of the propagation of light pulse in 
these active mediums have been theoretically discussed fl2l 
and experimentally demonstrated fl~3l l4l . In these experi- 
ments, the anomalous dispersion is usually realized by using 
a double-pumped gain medium in which the anomalous dis¬ 
persion lies in between the two gain peaks. As discussed by 
Kuzmich et.al ca, the gain medium is subject to quantum 
noise that accompanies the amplification process. In addi¬ 
tion, the gain medium could cause lasing when place inside 
a resonant cavity. To investigate how the quantum noise and 
associated gain influence the detector sensitivity and dynam¬ 
ics, we develop an input-output formalism for the optical field 
propagating through the gain medium. Using this formalism, 
we make a detailed analysis of the quantum shot noise lim¬ 
ited sensitivity for a typical gravitational wave detector con¬ 
figuration implementing the white light cavity idea, as shown 
in FIG. [I] Specifically, we consider: (i) the requirement for 
canceling the propagation phase shift; (ii) the optical stability 
of the interferometer system with the gain medium; (iii) the 
noise associated with the amplification process. Taking these 
factors into account, we find that the integrated shot noise lim¬ 
ited sensitivity is still limited by Eq.Q when the gain medium 
itself is stable (not lasing) and stationary. 


II. A BRIEF SUMMARY 


Before presenting the detailed analysis, we briefly summa¬ 
rize our main results in this section. The susceptibility of the 
double-pumped gain medium ^(Q) that we derive is given by 
(the same as in Refs. [12/14] but with slightly different nota- 




2 



(b) 



control fiold 




UJ h 


| 2 > 




signal field 
U p 



FIG. 1: (a) the typical dual recycled interferometer configuration 
for an advanced gravitational wave detector, with an atomic gain 
medium (blue block) embedded inside the SRC to compensate the 
phase delay of the arm cavity. An internal SRM (iSRM) with the 
same transmissivity as ITM is introduced to make impedance match¬ 
ing so that effectively we can view the compound mirror (consists 
of ITM and iSRM) as transparent to the sideband field fl6l : (b) the 
energy levels of the gain medium atoms. Two far-detuned strong 
control laser with frequency co a and co^ couple the energy levels |3) 
and 11). The signal field interacts with |2) and |3). 




FIG. 2: Phase angle and amplitude gain of the sideband field propa¬ 
gating through the atomic gain medium as functions of the normal¬ 
ized (by Aq) sideband frequency. The top figure shows the negative 
dispersion of the atomic gain medium. The white light cavity band¬ 
width is the linear region between — Aq and Aq. The bottom figure 
shows that the gain is negligibly small, except when Q ~ Aq. In 
these frequency regions, the gain is high and need to be considered 
in the design for preventing the possible instability (See Section III 
for detailed analysis). 


tions): 

y(Q) _ _ 2ir °P t _|_ 2ir °P t _ 
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( 2 ) 

where Ao is one half of the frequency difference between two 
control fields, and Q is the sideband frequency of the probe 
field with the carrier frequency G) p . The damping rate 712 is 
the effective atomic transition rate from state |2) to 11), while 
r opt , which depends on pumping power of the control fields, 
is the transition rate between |1) to |2) mediated by a virtual 
excitation of | 3 ). In terms of the ingoing and outgoing field 
d\ n , a oui are related by (temporarily ignoring additional noise 
term that will be mentioned later): 

a out («) = [l + ^(fl)/2]fl in (fl). (3) 

In obtaining x, we use the following approximation: 

Aq + (712 - r opt ) 2 > r 2 pt / 4 , ( 4 ) 

(For details, see Section III and Appendix B). Violation of 
this condition means the dynamics of the gain medium is 
non-stationary, as the two pumping fields start to interfere 
with each other. This approximation is also equivalent to 
the weak-coupling approximation which allows us to express 
the input-output relation Eq.Q for an unidirectional sideband 
field passing through the gain medium as: 

flout (Q) ~ e , * r(n)/2 e _ * i(a)/2 d in ( 0 ). ( 5 ) 

Here £ r (Q) and are the real and imaginary part of the 

susceptibility #(Q) of the medium, which describe, respec¬ 
tively, the phase accumulation and the amplitude change of 
the sideband field after passing through the medium as shown 

in FIG. H 


In order to compensate the round-trip propagation phase in¬ 
side the arm cavity thereby broadening the bandwidth of the 
optical cavity, the susceptibility should satisfy dx{Q)/d£l ~ 
—2L arm /c (negative dispersion), which leads to: 

r 0 pt[(ri2-r opt ) 2 -Ag] ^ E arm 

[(ri2-r opt ) 2 +A 2 ] 2 V' K) 

Once we fixed the parameter r opt and 712, we will have a pair 
of roots for Aq. For the positiveness of these Aq, the follow¬ 
ing condition has to be satisfied (see Appendix B for detailed 
derivation): 

(712 — ropt) < r opt c/( 8 L arm ). ( 7 ) 

Under these two conditions in Eq.([ 6 ]) and Eq.Q, we ex¬ 
plore the relevant parameter regime for studying the dynam¬ 
ical behavior of the gain medium. Firstly, as we analyze in 
detail in Section.III and IV(A), the system has two different 
types of instability (lasing) 1 ) if 712 < r opt , there will be a pop¬ 
ulation inversion between level | 1 ) and | 2 ), the gain medium 
starts lasing by itself, which we name as “atomic instability”; 
2 ) if the photon loss rate for each round trip inside the cavity is 
less than the photon increasing rate through the amplification 
by the gain medium, the cavity-medium system starts lasing, 
which we name as “optical instability”. In Fig[ 3 j we plot 
the phase diagram for the stability of the system. This figure 
gives a constraint on the possible parameter region for 712 and 
r op t of the atomic gain medium (with fixed SRM reflectivity 
r s ), if lasing were to be avoided. Notice that we choose the 
re-scaled parameter ( T] = r op t/7i2,^ = 8(712 - r 0 pt) 2 L/c7i 2 ) 
instead of (yi2,r opt ) and survey them within 0 < < 1. 

These new parameters help us exclude the atomic instability 
region (77 > 1) and the region where the phase-cancelation 
condition is unsatisfied (§ > 1). 
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FIG. 3: Stability region of the full interferometer scheme with 
double-pumped gain medium (optical stability only). The SRM 
power reflectivity rl = 0.5,0.8,0.9 are chosen from the top panel 
to the bottom, while we survey the parameter region for r opt ,yi 2 . 
The horizontal and vertical axis are rf = r opt /yi 2 and £= 8(712 — 
r op t) 2 T/cyi 2 . We survey 77 ,£ between 0 and 1 so that atomic insta¬ 
bility is excluded and the phase cancellation condition can be satis¬ 
fied. For each r s , the left panel and right panel correspond to two 
roots of Aq in Eq.{ 6 ), respectively. The purple region is the only 
stable region. In the blue region (“optical instability region”), the 
atomic medium is stable by itself but the dynamics of the full inter¬ 
ferometer system is unstable (see Section IV for details). The red 
region corresponds to the situation when the system becomes non¬ 
stationary, i.e. breaking down of condition in Eq.([4]). With the in¬ 
creasing of the SRM reflectivity, the stable region shrinks due to the 
enhancement of the optical instability effect 


Secondly, as implied by the above input-output relation, the 
gain medium is a parametric amplifier. Therefore, as first dis¬ 
cussed by Caves ED, there must be an additional noise term 
on the right hand side of Eq. 0 for keeping the commutation 
relation for <2 0Ut to be [a ou t(t), ^out(0] = 8(t — t'). This ad¬ 
ditional noise is due to the quantum fluctuation that causes 
spontaneous transition between |1) and |2). and degrades the 
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FIG. 4: Integrated shot noise limited sensitivity improvement factor 
(defined in Eq.(|45])) of the full interferometer scheme with double- 
pumped gain medium, taking into account of the effect of additional 
noise. The specification for the parameters is identical to the one for 
producing FIG. [3] The left panel and right panel correspond to the 
two roots of Eqjj 6 |. The dashed line is the boundary of the stable 
region shown in FIG. [3] It is clear from this figure that there is no 
parameter region where the integrated shot noise limited sensitivity 
improvement factor is larger than 1 , when the double-pumped gain 
medium itself is stable and stationary. 

signal to noise ratio. From the Hamiltonian, we can derive 
the noise terms from Heisenberg equations of motion. Their 
effect on the integrated shot noise limited sensitivity improve¬ 
ment factor (defined in Eq.(| 45 ])) is given in FIG. [ 4 ] (with tun¬ 
able parameters of atomic system and fixed SRM reflectivity). 

From these two figures, it is clear that 1 ) the stability con¬ 
dition and the phase cancelation condition put a strong con¬ 
straint on the possible parameter region; 2 ) There is no pa¬ 
rameter region where the shot-noise limited sensitivity is im¬ 
proved. This indicates that placing a stable double-pumped 
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gain medium with anomalous dispersion inside the SRC can 
not broaden the detection bandwidth while increasing the 
shot-noise limited sensitivity. Therefore, one shall explore 
other types of gain mediums or different parameter regimes 
for realizing the white light cavity. 


III. INPUT-OUTPUT RELATION OF DOUBLE GAIN 
ATOMIC MEDIUM 

After summarizing the main results, we now start a detailed 
discussion by first developing an input-output formalism for 
light propagating through the double-pumped gain medium 
in the Heisenberg picture. As we have briefly mentioned 
in the Introduction , our gain medium consists of three-level 
atoms schematically shown in FIG. [T] with two red (blue)- 
detuned (with respect to frequency difference between |3) and 
11)) control lasers. The polarizations of the control and probe 
fields are orthogonal to each other and only sensitive to the 
atomic transitions between |1) «->> |3) and |2) |3), respec¬ 

tively. In modeling the gain medium, we treat the atoms as 
non-interacting distinguishable particles. Nevertheless, all the 
atoms have the same energy level structures. In this section, 
we first derive the atomic dynamics for a single three-level 
atom, then extend the result to many-atoms case under the 
approximation that the length of the gain medium is much 
smaller than the spatial scale of the optical sideband field 
2 nc/£l where Q is the gravitational wave frequency. 


A. Single-atom dynamics 

The above physical modeling leads to the following system 
Hamiltonian for a single atom 

H = H atom +H f + H int +H r ( 8 ) 

The //atom is the free Hamiltonian for a three level atoms in 
the form of: 


F/atom — 


O aa ■ 

a= 1,2,3 


(9) 


where (Q a is the Bohr frequency of energy level a and 6 aa is 
the atomic population operator. 

The Hf is the free Hamiltonian for the sideband probe fields 
propagating in the main GW detector. Since we are only in¬ 
terested about optical modes that are around the central fre¬ 
quency of the probe field co p , these modes have frequency 
co p ± £1 where Q denotes the frequency band that we are fo¬ 
cusing on. Then we have the Hamiltonian Hf as: 


H f = he 
« he 



- k')a f _ kp+k ,a_ kp+kl 

— k )& r _ k p +k ,a_ kp+k i , 


( 10 ) 


where we have assumed that the field propagate uni- 
directionally (from right to left) and k p = (O p /c,Ak = £l/c. 


We have also used the narrow band approximation A k/k p = 
kl/co p <C 1 so that we can extend the integral range to 
[— oo ? oo]. By defining the optical creation/annihilation oper¬ 
ators in x—domain as: d x — f^dk'a_ kp + k re~ lkx , the above 
Hamiltonian can be written as (for details, see [20, 21 ] or Ap¬ 
pendix A)\ 


Hf 



dd\ 

~d^ C 


- da x 
dx 


dx. 


(ID 


Notice that the a x is the slowly varying amplitude operator 
(both spatially and temporary) with respect to ? 

defined as: E p (t,x) = a x (t)e~ mpX l c ~ lGipt +h.c. The probe field 
propagates unidirectionally (right to left), encounters and in¬ 
teracts with the atom at position xq. This interaction is given 
by a Jaynes-Cumming type of Hamiltonian under the rotating 
wave approximation [18]: 

Hint = - 2 ^dl Q e l(0pt 023 - -IIysE*Gi$ + h.c, (12) 

where the first term describes the atomic transition between 
12) and |3) under the driving of probe fields with transition 
operator 623 and the second term describes the atomic transi¬ 
tion between |1) and |3) under the pumping of control fields 
with the transition operator 613. The E c = E a e l(0at + E^e 1 ^ 
describes the two classical amplitude of control fields with 
frequency CO a f and the ii mn (m,n = 1,2,3) are the dipole mo¬ 
ments of the atom. The atom transition operators satisfy the 
algebra: 


&mn&kl — Gmldyik > ( &mn )^ — Onm- ( 13 ) 

The coupling between an atom with other bath sources at 
position xo is introduced phenomenologically by: 

H y =ih^/ 2 Y\ 2 n\ 2 e lG>llt 6\2 - ^i3ale l(0ot O 13 +h.c (14) 

where fin is the noise operator which couples to the atomic 
transition operator between |1) and |2). The a c is the quan¬ 
tum fluctuation associated with the control field. In our 1- 
D model, the quantum fluctuation associated with the probe 
field has been included in the a Xo field of Eq.([T2]). The noise 
bathes model can be attributed to mutual collision of atoms or 
the stimulation of the external electromagnetic vacuum bath. 
Here to study the minimal additional noise, we consider an 
effectively zero-temperature external bath. 

With the Hamiltonian, we can now analyze the dynamics of 
the gain medium. The Heisenberg equation of motion for the 
probe field derived from Eq.([TT]) can be written as: 

d ~t~ cd ~t = (15) 

which reflects the fact that the probe field propagates from 
right to left (unidirectional). 

The Heisenberg equations of motion for the atomic transi¬ 
tion operators of a single atom in this gain medium are given 
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by: 


£ 0 C) dynamics can be written as: 


613 + (/'©31 + 713)613 = iy/ 2 Yi 3 (ou ~ <733 )a cm e i(0ot 

+ ^13(611 - d33)E c + -i*nana Xo e- im i* ( 16 a) 

623 + (ia >3 2 + 723)623 = ^13621 (E c +a c e~ i(0ot ) 

+ *0/^3(622 — 633)a in e _ '“ pf ( 16 b) 

612 + (ia>n + 712)612 = -\/27u(6ii - 6 2 2)«i2ine“ ,a)21? 

^23613^0- ^13632(£ e +a c e“ !<a3 0 ( 16 c) 

Notice that the fti2in> <^in and a c i n are the incoming noise fields, 
whose relations with the n\2,a,a c of Eq.(p 4 |) are given in the 
way ofEq.@£l]). The 713 = Mi3/^?23 = M23/8 can de- 
rived from Eq.(| 12 |). Besides, the condition that the majority of 
atoms are initially prepared at | 1 ) is set as an assumption. In 
a real experiment, this population preparation can be achieved 
through various methods such as introducing an additional op¬ 
tical pumping field fT4l . 

It is clear that the above equations of motion are generally 
nonlinear. However the system dynamics can be simplified 
by exploring the linear regime where the scheme is proposed 
to operate. The simplification can be done using perturba¬ 
tive method by noticing that 1) the control fields have large 
detuning with respect to CO31 and therefore the population of 
atoms on 13) is still small compared to that on 1 1); 2) the tran¬ 
sition between 1 1 ) — |3) is much stronger than the transition 
between 1 1) — |2) and | 2 ) — |3) since it is induced by strong 
control beams; 3) The probe field is very weak compared to 
the control field since it is around the quantum level. 

There are three dimensionless expansion parameters in this 
system of equations of motion: £ ~ fi mn E c /A 0, a ~ H mn a /Ao 
and a < £ < 1 (notice that the denominator could also be 
other frequency scales such as 0)31 — (D a )b , we choose the 
smallest one here for briefness). Writing the 613, 623, 612 
in the rotating frame of (Do = ((D a + (D b )/ 2 , (D p and (Do — (D p 
respectively, the leading order £) of 613 dynamics can be 
derived as: 


613 —i(fflo-©3 i + *713)613 = ^13611 ( E a e iA °‘+ E b e iA °‘), 

( 17 ) 

in which we can approximate the collective population op¬ 
erator on 11) as an ~ 1 and a mn is the mean value of o l mn 
averaged over all the atoms. The solution of Eq.fTT]) is given 
by: 


623 -i((D p — (D32 + *723)023 » 2^ 13< ^ 21j ® c ’ (l^a) 

a „ i _ i „ ~ 

012 + 712012 = 2 A^23^x 0 013 - 2^13032^ 

- in, ( 19 b) 

where E c = E a e~ lAot + E b e lAot is the pumping field strength 
in the rotating frame of Ctfo and we have used the fact that 
(Do = (D p + U)2i (See FIG. [I]). We also make use of the fact 
that O', 7 , = 1 . In deriving Eq.(p~ 9 ]), we also assume that system 
parameters satisfy: 723 /{cd p — CO32) <C £ 2 0 C. 

For the probe field, we can integrate Eq. © around xo and 
obtain: 


1 

- 4 0 + + 4 0 - = 2 ^23 <723, (20) 

in which the a Xo + and a XQ - are the incoming and outgoing 
sidebands fields (respectively) defined in the vicinity of the 
interaction point vo (in the following, we will use a m and d out 
to represent them, respectively). The probe field at xo is con¬ 
nected with these vicinity fields through the junction condi¬ 
tion: 


d*0 — 2 (^ n ^out) • (21) 

The dynamics of <732 can be obtained by solving Eq.(p~ 9 ])(a), 
we have: 


032 = 


M13 K A 

—- 012 , 

2 (D32~CD p 


( 22 ) 


In deriving the above equation, we assume that 723 <C 0)32 — 
(D p and make use of the fact that O12 is a slowly varying am¬ 
plitude thereby 632 ~ 0 in Eq.p 9 | ). Substituting the solution 
Eq.(| 22 |) and Eq.(| 2 T]) into Eq.( 19 |)(b) and E q ([ 20 ] ), we can adia- 
batically eliminate the 623 so that the Eq. ( 1911 (b) and Eq.pO]) 
form a closed equation set: 


4 ut = ain - KV^Yo pto e ‘ A °' + \/ 2 7 optb e“ !A °')612 ( 23 a) 

612 + 712612 = i (v /2 7 o P tbe“' Ao< + ^roptae^) aj n 
+ ( 7 opt + 72 Ao) 6 i 2 - /© op t6i2 - *j 2 Ynh n m- ( 23 b) 


This equation set describes the coupling between the compos¬ 
ite system (consisting of the atom and the pumping fields) to 
the probe field. 

The second term on the right-hand-side of Eq.(| 23 ]) is the 
sum of an anti-damping term y op t 012- 


013 


1 M13 E a e-“* 

2 0)31 -(Da 


2 (D3l-(D b 


(18) 


Here, we have neglected the 713 which is assumed to be much 
smaller than the detuning: 713 /{(D a j, — 0)31) £. In the same 

rotation frame, the leading order of the 623 (^ £ 2 0 C) and 612 


7opt — Yopta + Yoptb — 


M 2 2 3 Mi 3 r 
32 


\E a \ 2 | E ,' 2 


A l 




A* 


and a high-oscillating term 72 a o 012- 

^23^13 


72A 0 = 


2 " 2 v E a E p e 2iAot E b E*e~ 2 iA o* 


32 


A1 


A1 


( 24 ) 


( 25 ) 
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These formulae are derived under the approximation Ao 
A p ^b (thereby A p « A a « A^). These are good approximations 
to the situation in the proposed experiments EBnaiBi. 

In the symmetric pumping case when E a = Eb = Eq: 

7opt W M23Ml 2 3^o/( 16A «)- (26) 


When y opt is larger than 712 , we have the population inversion 
between energy level | 1 ) and | 2 ), i.e., the atomic instability, 
mentioned earlier. 

Solving the Eq.(|23])(a)(b) in the frequency domain, we can 
obtain the input-output relation for the probe field: 


flout(£i) =^(^)flin(^) + <^+(^)«{ 2in (Ao - &) 

+ ^_(O)4 in (-A 0 -£>), 
with ^(Cl) and ^(£1) given by: 


(27) 


= 1 - 


7opt 


7opt 


i(£l + Ao) — 712 + 7opt ~ Ao) — 712 + 7opt 

(28a) 

\/2yi27opt 


^±(Q) = 


±/Aq - i£l + 712 - 7opt 


(28b) 


Notice that 1) <A±(—£l) = 2) Here and after, for 

simplicity, we will only consider the symmetric pumping case 
where E a = E\ because the non-symmetric pumping will only 
induce an additional rotation int he quadrature plane, which 
does not affect our main results. 3) In obtaining the required 
susceptibility Eq.Q and the corresponding input-output re¬ 
lation, we neglect the above oscillating terms at 2 Ao which 
is contributed by the driving from the beating of two control 
fields. This approximation leads to the condition in Eq.([ 6 ) 
(Details are discussed in Appendix B); 2) The tiny Stark fre¬ 
quency shift G) opt = (i \3 |£ c (f)| 2 /(4A,) Ao on the right hand 
side of Eq.([23]) has been neglect here. 

The above input-output relation describes a phase- 
insensitive parametric amplification process. Therefore, there 
is an additional noise given by the fij 2in terms in Eq.(|27|). This 
noise comes from the stochastic dynamics of 612 driven by the 
last term on the right hand side of Eq.(|23]). The usual method 
for introducing the additional noise for parametric amplifier 
is using the argument given by Caves ED, which base on the 
principle that a out field should satisfy Bosonic commutation 
relation. However, Caves’s method can not be directly applied 
to our system since the additional noise has two different fre¬ 
quency channels. Solving the dynamics from the full Hamil¬ 
tonian of the system can allows us to pin down the source of 
the additional noise and give the correct formula for the noise 
contribution. 


B. Extension to many-atoms case 

In the above subsection, we discussed the input-output rela¬ 
tion for the probe field interacting with a single atom. In this 
subsection, we extend the above results to the many-atoms 
case. 


Since the size of the gain medium (centimeter scale) is 
much smaller than the spatial scale of the slowly varying 
amplitude of the probe field (kilometer scale), therefore the 
slowly changing amplitude of the probe field interacts with 
all the atoms together. In this case, the anti-damping rate 
will be enhanced by a factor of N where N is the number of 
the atoms [22 -24I and the ^ coefficients in the above input- 
output relation can be written as: 


= 1 - 


r op t 

f(Q +Ao) — 7i2 + r opt 


r op t 

i(£l — Ao) — 712 + r opt 
(29a) 


where r opt = Ny 0pt . 

The formulation of noise field in the input-output relation 
for many-atoms case has some subtleties, depends on the spe¬ 
cific modeling of the interaction between the noise field and 
the atoms. 

• Noise interacts with atoms locally —In this case, each 
atom is associated with its own noise bath. The noise term 
will be represented by: 


N 

Noise term = £^^ ± (0)«f 2 in (±Ao-n), (30) 

±7=1 


where 


oY±(G) 


y/27l27opt 

±/A 0 - i£l + 712 - r opt 


( 31 ) 


• Noise interacts with atoms collectively — In some cases, 
the noise is introduced through processes where the electro¬ 
magnetic field amplitude interacts with all the atoms collec¬ 
tively as what the slowly-varying amplitude of the probe field 
does. For example, the 712 is induced by applying an addi¬ 
tional pumping laser such as the experiment done in (TO). In 
this situation, the noise term will be represented by: 

Noise term = (^)fij 2 in ( c ° 21 ± Ao — Q), (32) 

± 


with 


JE£(Cl) = 


V 2 7L2r 0 pt 

zh'Ao — + 712 — r opt 


(33) 


Notice that the 712 here (also accordingly in ^(Q)) should 
be understood as N times of the transition rate from |2) to 11) 
for one single atom which is proportional to the intensity of 
the additional pumping laser the experiment in Go). 

Notice that 1): The input-output relations based on both of 
these noise models satisfies the Bosonic commutation relation 
[a 0 Ut (fl),< 2 0 Ut (fl / )] = 5(Q — Q!) under the weak coupling ap¬ 
proximation (See Appendix B). For the single-pumping case 
where Aq = 0, the Bosonic commutation relation will be ex¬ 
actly satisfied. 2): More importantly, as we shall see later, the 
subtleties of the noise model do not affect the sensitivity of 
the gravitational wave detector. 
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C. Some physical discussion 

After deriving the system dynamics and input-output rela¬ 
tion, we are going to give some intuitive discussion of the sys¬ 
tem dynamics and the additional noise. 

Firstly, the “anti-damping” dynamics of 612 can be under¬ 
stood in the following way: a small amount of atoms initially 
populated on | 1 ) can be pumped to |3) by the detuned con¬ 
trol fields, then jump to | 2 ) due to their interactions with the 
probe field. During this indirect transition between |1) and 
1 2 ) mediated by |3), the population of atoms on \ 2 ) will in¬ 
crease indefinitely if the decay rate from | 2 ) to | 1 ) is not suf¬ 
ficiently large-a “population inversion process”. Physically, 
this process could cause lasing (“atomic instability” in Sec¬ 
tion II) and our approximation will fail as the population on 
12) becomes larger than the population on 11). One may argue 
that this instability will not happen in real experiments with 
the atom population being prepared using additional pumping 
fields. However, the thermal collision relaxation rate can be 
tuned to be small if we decrease the gas temperature, increase 
the pumping beam waist and fill in the “buffer” gas [J_2, 19]. 
In this case a small transition rate contributed by the optical 
pumping beam could be sufficient for the population prepara¬ 
tion. Therefore, in principle the lasing can still happen as long 
the control fields are strong enough and r op t > 712 . 

Secondly, for the additional noise, the stochastic dynam¬ 
ics driven by the h\ 2 m can be attributable to 1 ) the collision 
of atoms due to the Van-der-Waals mutual interaction or the 
thermal collsion m, 2 ) the transition between | 2 ) — 11 ) in¬ 
duced by environmental black-body radiation, 3) the contribu¬ 
tion of the quantum noise associated with the additional opti¬ 
cal pumping process as in ifTDl . In Eq.@7) and Eq.@5)-g2), 
the a out field contains terms related to the additional noise h 
in the way that the stochastic fluctuations of the population on 
| 1 ) and 12 ) will cause fluctuations of the transition between 
1 2 ) and |3), since 62 3 is slaved by 612 . 

In this section, we have derived the input-output relation 
for the sideband probe field propagating through the double 
gain medium from the full Hamiltonian. We also discussed 
the opto-atom dynamics and origin of the additional noise. In 
the next section, we will apply these results to the interfer¬ 
ometer configuration shown in FIG. [T] and analyze its strain 
sensitivity. 


IV. INTERFEROMETER WITH GAIN MEDIUM 

The propagation of sideband fields inside the interferom¬ 
eter shown in FIG. [T] can be schematically described by the 
flow chart shown in FIG. [5] Here, we only study the differen¬ 
tial mode of this interferometer which carries the gravitational 
wave signal and can be mapped into a signal cavity contain¬ 
ing a gain medium [25]. In this scheme, an internal signal 
recycling mirror is used to effectively remove the frequency 
response of the arm cavities so that the sideband fields see 
an input-output relation for the amplitude and phase quadra¬ 
ture of sideband field in the following form (we ignore the 
optomechanical back-action term by assuming the infinitely 


GW 



FIG. 5: A flow chart showing the propagation of fields in the full 
system. The optical stability of the system is determined by the part 
inside the blue dashed box, whose gain function is given by Eq.(??) 
and ( [44] ). The test mass (end mirror) m is driven by gravitational 
waves while the double gain medium is affected by the additional 
noise h. 

heavy test masses): 

e(Q) = M 0 (Q)d(Q)-hD(Q)h(Q), (34) 

or more explicitly, we can expand the vectors d,e,D and ma¬ 
trix M 0 


' w 


'1 o' 

’rfi(Q) " 

+ e iaz 

0 

. h(&) _ 

— t 

0 1 

. d 2 (a) _ 

_ V2j?_ 


(35) 

with = P c (OoL 2 /(he 2 ). 

(A) In case of local noise model, the input-output relation 
for phase and amplitude quadrature of light field propagates 
through the gain medium is (Ref. Eq.([29])-([32|)): 

d(Q) =M(Q)c(Q)+N+(Q)+N-(Q). (36) 

where N± represent the additional noise terms, the forms of 
which depends on the specific noise modeling as we will show 
later. 

The combined effect of the gain medium and the main in¬ 
terferometer is described by M tot (Q) = Af (Q)Mo(Q) with the 
noise terms N±(£l)Mo(£l). Then the final input-output rela¬ 
tion for the sideband field is given by: 

j(G) = M k (£l).k(Q.)+t s e iaz M c (Q.).D(Q.)h(£l) 

+ t s e 2iaz M c (Q.).(N + (C1) +N-(£l)), *"' 

in which Mfc(£2) = — r s I + ^M c (Q)M tot (Q) and M C (Q) = 
[I — r s M to t(^)] _1 with r s ,t s the amplitude reflectivity and 
transmissivity of the signal recycling mirror. The k(Q) and 
j(Q) are the input and output field of the entire configura¬ 
tion as shown in Fig [I] The first term in Eq.([T7]) is the quan¬ 
tum noise contributed by the vacuum injection outside of the 
SRM, the second term is the contribution of additional noise 
introduced by the gain medium and the last term is the sig¬ 
nal term. The SRM feeds the optical field back into the gain 
medium and the main interferometer, which is described by 
dashed blue box in the flow chart FIG. [5] This feedback pro¬ 
cess will bring in another potential lasing (the “optical insta¬ 
bility” mentioned in Section II) even if 712 > y op t, discussed 
in the following subsection A. 
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For the additional noise terms N±(£l) in the above rela¬ 
tions, suppose the additional noise interacts with atoms lo¬ 
cally, they are given by 

N 

N+(Cl)+N-(Cl) = £[N ± {£l).n j± {£!)}. (38) 

M 


The are the vectors (/rf^-n^), with ^ is the am¬ 
plitude (phase) quadrature of the noise field n J l2in with respect 
to the central frequency ±Ao. The noise matrix N± can be 
derived using the sideband-quadrature transfer matrix M^ v de¬ 
fined as: 

= !)• <39) 

which leads to: 

N ± (f!) = M„. ( ^ ^° (n) ) ,!«„> (40) 

in obtaining the above formula, we have made use of the rela¬ 
tion: 

If the noise interacts with the atoms collectively, the noise 
term is given by: 


where formally we have: 


•M- 1 n ± (n) 

(41) 


<A±(Q.) = VN 


\/27i2y 0 pt 

±/Aq — + yi2 — r opt 


(42) 


where we have used the fact that r opt = A/y op t while as before, 
the 712 here should be understood as N times of the |2) —>> 
|1) transition rate for a single atom. When we calculate the 
sensitivity, the y/N in the above equation will give us the same 
N factor in the following Eq.(|43]) while the rest part of the 
above equation has the same form as JV±j. Therefore if the 
the noise is modeled as interacting with the atoms collectively, 
the strain sensitivity formula will have the same form except 
that the 712 are different. However, since we will scan over all 
possible values of 712 and r opt in the sensitivity calculation, 
this difference of 712 will not affect our final result. 

Finally, the shot-noise limited strain sensitivity of the inter¬ 
ferometer derived from Eq.(|43]), which is given by: 




kiV/,.M c (Q).D(n)| 2 

v*M c (J2).N J - ± (fl).Nt t (J2)Mt(J2)vJ 
¥ |v/,.M e (Q).D(n)| 2 

(43) 


in which the V/* = (sin £,, cos £,) is the homodyne readout vec¬ 
tor with £ is the homodyne angle. In our calculation, we will 
choose § = 0 which means we only measure the phase quadra¬ 
ture 72 (£2) of the output field. 

We will discuss the numerical result of this shot-noise lim¬ 
ited strain sensitivity in the following subsection B. 


A. Stability Criteria 


Besides the atomic instability due to “population-inversion 
process” when 712 < r opt , it is important to notice that the 
dynamics of the interferometer with the gain medium may still 
be unstable (start lasing) even if 712 > r opt . 

This instability is related to the feedback process discussed 
below Eq.([37]) due to the reflection of the SRM. Intuitively, 
when the reflectivity of SRM r s becomes high (or equivalently, 
t s decreases), the photon loss rate through the transmission for 
each round trip can be less than the photon increasing rate 
through amplification by the gain medium, corresponding to 
the “optical instability”. The criterion of this instability is 
determined by the analytical behavior of the close-loop trans¬ 
fer matrix M c . In our configuration, this close-loop transfer 
matrix is diagonal so that it can be simplified as a close-loop 
transfer function G C (Q) (M c = G C (Q)I): 


G C (G) 


1 

1 -r s G 0 (Q) 


(44) 


with G 0 (Q) = e 2lQr M(Q). The stability of the full system 
is determined by the poles of the denominator, which can be 
obtained by solving the equation 1 — r s G 0 (Q) =0. 

However, the time-elapsed factor e 2lQ,T in the gain func¬ 
tion makes it difficult to find the root of the above mentioned 
equation. The Nyquist criteria provides us another way to 
understand the stability through the analytical behavior of 
G 0 (Q) [27 ](See Appendix C) instead of G C (Q). Specifically 
in our system, the Nyquist criteria can be stated in a way that 
the Nyquist contour of r s G 0 {0 .) should not encircle the point 
(1,0) in the (Re[r s G 0 ],Im[r s G 0 ]) plane at all. This criteria is 
equivalent to the lasing condition that the round-trip gain is 
smaller than one when the phase is the integer number of lit. 
For illustration purpose, several examples of the Nyquist con¬ 
tour of r s G 0 {£l) are shown in FIG. [b] given typical parameters 
of 712 and Ao. This plot demonstrates that increasing the SRM 
reflectivity can lead to system instability. We further search 
the parameter region 0 < 77 , ^ < 1 and give the plot on Fig[3j 
from which we can see that the stability criteria imposes a 
very strong constrain on the possible parameter region. Only 
the interferometer system with parameters in the stable region 
is useful. 


B. Integrated shot noise limited sensitivity improvement factor 

For quantitatively describing the improvement of the inte¬ 
grated shot noise limited sensitivity, we define a quantity: the 
integrated shot noise limited sensitivity improvement factor 
(iSNS improvement factor) p r in the following way: 



where S% h is the shot noise limited gravitational wave strain 
sensitivity of the laser interferometer with double-pumped 
gain medium given by Eq.(|43]), while the denominator is that 
of the conventional interferometers given in Eq. ([T}. In case of 
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- 1.0 - 0.5 0 0.5 1.0 

Re[r s G 0 ] 


FIG. 6 : Nyquist contours of the r s G 0 (£l) for the full system with 
fixed parameters of the gain medium 17 = 0.1,£ = 0.4 while vary¬ 
ing the SRM amplitude reflectivity r s . The dashed (magenta), solid 
(red), dotdashed (black), dotted (blue) curves are the Nyquist con¬ 
tours when r% = 0.9,0.8,0.7 (unstable cases), and 0.5 (stable cases), 
respectively. The upper-part zooms-in the full contour (lower-part) 
near (1,0) (the red spot). It is clear here that when the SRM reflectiv¬ 
ity increases, the instability develops and the stable region in FIG. [3] 
shrinks. 


p r > 1, the system with double gain media will improves the 
signal to noise ratio by breaking the trade-off between the de¬ 
tection bandwidth and peak sensitivity. However, we need to 
be cautious about the stability of the system at the same time. 

We can calculate the strain sensitivity and hence the iSNS 
improvement factor. By fixing the SRM power reflectivity to 
be rj = 0.8, we calculate the iSNS improvement factor by 
searching the parameter region for (tj,£) within the range 
[0,1] constrained by the phase cancelation condition. For il¬ 
lustration purpose, we first calculate the p r by ignoring the 
effect of additional noise introduced by the atom system and 
give the plot in FIG. [7] This figure clearly shows that there 
could be some parameter regions where the system is stable 
and p > 1. 

However, taking into account of the additional noise, the re¬ 
sults dramatically changed as we can see from FIG. [4] It turns 
out that there is no region where p > 1 and the p — contours 



FIG. 7: Integrated shot noise limited sensitivity improvement factor 
(defined in Eq.(|45|)) of the full interferometer scheme with double 
gain medium, without the effect of additional noise. The specifica¬ 
tion for the parameters is identical to the one for producing FIG. [3] 
and FIG. [4] The left panel and right panel correspond to the larger 
and smaller roots of Eq. & respectively. The dashed line is the 
boundary of the stable region shown in FIG. [3] In this figure, when 
the detuning takes the larger solution, there are some regions where 
p > 1 and the system is stable at the same time. However, as we 
can see from FIG. [4] these regions will disappear when we take into 
account of the effect of the additional noise. 


are significantly distorted due to the additional noise. Accord¬ 
ing to our numerical test, this conclusion does not change with 
the variation of the SRM reflectivity. 


V. CONCLUSION 

In this paper, we applied the input-output formalism devel¬ 
oped from the Hamiltonian of light-atom interaction to study 
the quantum noise of white light cavity using double gain 
medium. We find that not only does the additional noise as¬ 
sociated with the parametric amplification process affects the 
system, but the requirement for the system stability also in¬ 
troduces an additional issue to take into account for its im¬ 
plementation. We conclude that the net sensitivity can not be 
enhanced by using the anomalous dispersive behavior of the 
stable double gain medium when the system is stable. For 
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further study, we will consider the situation that the system is 
unstable but being controlled by a external feedback loop in 
an accompanied paper. 
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Appendix A: Slowly-varying amplitude Hamiltonian of the 
electromagnetic field 


In the main text, the Hamiltonian of the electromagnetic 
field is given in Eq.fTT]). Unlike the usual free field Hamil¬ 
tonian written in the k— space, this Hamiltonian is written in 
the v—space and the a x is the slowly varying amplitude of the 
optical field. In this appendix, we give an derivation of this 
form of Hamiltonian. 

In the k— space, the free-field Hamiltonian for an unidirec¬ 
tional propagating field can be written as: 

fAk 

Hf = hcJ dk'(k p -k')a\ kp+k ,a_ kp+k i , (Al) 

Since we are interested in the Ak <C ko, we can approximate 
the above formula to be: 

H f « he J dk'(k 0 - fc , )a! jto+jt /a-* 0 +*' • (A2) 


This approximation is called narrow band approximation. 

Then we can define the optical field operator in the v—space 
by Fourier transformation: 


a x m / d _£ 0+jt /e lk ' x dk’ 

J —oo 


(A3) 


Substituting the above definition into the Eq.( |A2| ), we obtain: 
ihc 


Hf = 


= heko J 


dxala x + ^ 


j: 


dx 


dal 

lh C 


c-fda x 

7 - 

* dx 


(A4) 


Further more, if we work in the rotating frame of cko: a x = 
d x e~ lck °' then the above Hamiltonian will be: 




djl 

dx' 


da x 

dx 


(A5) 


Notice that the d x satisfies the commutation relation: 
[a x , dj} = 8 (x — x r ) . The a x here is the spatially and temporary 
slowly varying amplitude of the electromagnetic field. This 
fact can be seen from the definition of the electric field under 
the narrow-band approximation: 

P + \x,t) = E^ + \x,t)e i(0 ° ,+ik ° x 

sa \j2nhckf) I k-Ai c dk'a_i C0 -k'e^ x ~ lck ' t = \J2ntick{)d x 

(A6) 
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Appendix B: Validity condition and Phase cancelation condition 


In the main text, we mentioned that %(Q) given by Eq.([2]) 
and the input-output relations Eq.([27]) and Eq.([29|)-([32|) are 
obtained by the neglecting the terms with frequency ~ 2 Ao. 
Here, we are going to give a derivation of this term to show 
how it affects the field input-output relation of the double gain 
medium. The order of magnitude of this effect determines the 
validity region of Eq.( [29a| ). Notice that for simplicity, we only 
discuss the case of symmetrical pumping with E a = = E c . 

We start from the single-atom dynamics. As we can see in 
the main text, the equation of motion for 612 after the adia¬ 
batic elimination of 623 can be written as: 

612 + (712 - 7opt)<7i2 =^23&lo n + 

z (El) 

- i2in- 


Here for simplicity, the pumping strength is assumed to be 
symmetric: 


<7i 3 (» ~ ~ l?>Ec cos(A 0 1). (B2) 

CO 31 -CO 0 


Assume that d \2 —» <7® + where of^ is the part which 
satisfies the Eq.( |Bl| ) without the 72 a 0 term and off is due to 
the effect of 72 a 0 term. Then using iteration method, off sat¬ 
isfies (temporarily neglecting the ft 12 term): 

<*12 + (712 - 7 opt )= fiAo • (B 3 ) 

The djj in the above equation can be solved as: 

/ e^ n ~ ro ^^~^ cos(Aot')dl n (t')dt', 

J -OO 

(B4) 

which leads to: 


72a ( A ( 2 ) 0) — 2 ToptVTtpt J d£lal n (£i) 


g—i(Ao+Q)t 


J{A 0 -Cl)t 


+ 


712 - 7opt - i( Ao + £2) 


7i2 — 7opt + i(Ao - £2) 
+ O(3A 0 ). 


(B5) 


The terms oscillating around 3Ao will be neglected in the fur¬ 
ther discussion. Substituting the above formula into Eq.( |B3| ), 
we obtain off d as: 


«T(«) = 


z’v%pt7opt 




2(7i2 - 7o P t - *£2) 

4n(£2-A 0 ) 


7i2 — 7opt — *(£2 + 2A 0 ) 


+ 


7i2 — 7o P t — *(£2 — 2 Aq) 


(B 6 ) 


Based on this equation, using the relation Eqs.(20) and (22), 
we have the input-output relation as: 


M(fl) =1 -hf-(Q) 


_7>pt_ 

—iQ — /Aq + 712 — 7opt 


+ /+(^) 


_7opt_ 

—iQ + /Aq + 712 — 7opt 


(B7) 


in which 


f±(Q) — 1 + - 


7opt 


2 712 - 7opt + ^^^o) ’ 


(B8) 


where y„ pt is given by Eq.([24|) of the main text. Then the va¬ 
lidity condition for our input-output relation is clear: 

|/±-l| 2 «l. (B9) 


In the balance pumping case and under the D-C assumption 
Q ~ 0, we can simplify the above validity condition to be: 

Aq + (712 - 7o P t) 2 » (7o P t) 2 /4. (BIO) 


In the case of many-atoms case, the t in the above valid¬ 
ity condition will be replaced by r opt . 

If the validity condition for the input-output relation was 
satisfied, the input-output relation Eq.([29|)-([32|) can be used 
in analyzing the system. Moreover, as we have mentioned in 
Section I, this condition is equivalent to the weak coupling 
approximation which allows us to approximate the phase gen¬ 
erated by the gain medium to be Re[#(Q)]/2. Besides, for 
achieving cancelation of the propagating phase inside the arm 
cavity in the weak coupling limit, another condition must be 
satisfied: 

r op t[(7i2 ~ r op t ) 2 — Aq] _ _ L arm 11 ) 

[( 7 l 2 -r opt ) 2 -A 2]2 - c ' 

This condition will reduce to r opt = AqL arm /c when |yi 2 — 
r opt| <C Aq. In our calculation, we have used the exact formula 
Eq.(BTTT). 

If we fixed the value of 712 and r opt , then the phase cancela¬ 
tion condition becomes a second order algebraic equation for 
Aq. Suppose this equation has two roots xi,X 2 , then we have: 


X 1 X 2 — (712 r^) [(712 r opt ) + r opt c/( 2 L arm )], 
x\ + X 2 = r opt c/( 2 L arm ) — 2 ( 712 _ r^) • 


Notice that in the above equations, X\X 2 is always positive, 
thereby x\ +X 2 can only be positive: 

(712 — r opt ) < r opt c/(4L arm ). (B13) 

On the other hand, Eq.( |Bll| ) must have real roots, which 
gives: 

(712 — r opt ) < r opt c/(8L arm ), (B14) 

which is a more stringent condition than Eq.( |B13| ). 

In summary, considering the validity of our input-output 
relation Eq.(29|)-([32|) and the requirement of the phase- 
cancelation condition, our parameters must satisfy Eq.( |B10| ) 
and Eq.( |B14| ). It is important to notice that there are always 
two Aq corresponding to a fixed set of ( 7 i 2 ,r opt ). In plotting 
the FIG.[4]and FIG.|7j we should take into account both roots. 


























12 


Appendix C: Nyquist criteria 

In Appendix C, we will give a brief introduction to Nyquist 
criteria ma which we used in understanding the stability con¬ 
dition of the full system in Section IV. 

The behavior of control systems is usually described by 
gain functions. For a control system with feedback pro¬ 
cess, the open-loop gain function G 0 (Q) is used to describe 
the information transfer ignoring the feedback process, while 
the closed loop gain function G C (Q) includes the effect of 
the feedback process. The relationship between G 0 (Q) and 
G C (Q) can be written as: 


G C (Q) = 


Go (£1) 

1 +H(Q)G 0 (Q)' 


(Cl) 


The H(Q) is the gain function for the feedback process itself, 
it is clear from FIG.[5]that in our system, it is just the reflection 
of the SRM: — r s . 

The stability of the system depends critically on the poles 
of the close-loop transfer function, that is, it depends on the 
poles of G 0 (Q) and also the zeros of 1 — r s G 0 (Q). However, 
computing the poles and zeros of these gain functions is gen¬ 
erally a difficult task when these functions are non-rational. 
The Nyquist stability criterion is a graphical technique for de¬ 
termining the stability of a control system, which is based on 
the following Lemma : Cauchy argument principle. 

The Cauchy argument principle starts from the Nyquist 
mapping, which maps the complex argument fl-plane to the 
complex F(Q) -plane. If we have a clockwise contour in Q- 
plane encircling a zero of F(£l), correspondingly, the con¬ 
tour also encircles the origin clock-wisely in the F(Q)— plane. 


However, if we have a clockwise contour in £2-plane encir¬ 
cling a pole , then the corresponding contour will encircles the 
infinity clock-wisely in the F(Q) —plane thereby encircling the 
origin in an anti-clockwise way. In general, if we have a con¬ 
tour in the F(Q) -plane encircling the origin N times clock- 
wisely, that means in the Q—plane, the number of zeros (Z) 
and the number of poles (P) satisfy: 

Z = N + P. (C2) 

This equality is the Cauchy argument principle. 

The transformation for quantity A(t) between the fre¬ 
quency domain and the time domain is defined as: A(t) = 
f^° OQ A(Q)e~ lQt . Therefore if A(Q) has poles in the upper- 
half plane, we will have instabilities for casual system (t > 
0). Now we choose the contour encircling the upper-half 
Q— plane as “Nyqusit contour”. If the system is stable, then 
the Z of 1 — r s G 0 (Q) (the denominator of close-loop gain 
function) inside the Nyquist contour should be zero. As a re¬ 
sult, the Cauchy argument principle becomes N = —P, which 
is the Nyquist criteria. 

In our system with G 0 (Q) = %(Q)e 2lQT , we have the the 
poles of r s G 0 (Q ) which are Q = ±Ao — i(yn — r opt ). Both 
of them fall outside of Nyquist contour because yn — r op t > 0 
for the requirement of the stability of the atomic gain medium 
itself. Then we can conclude that P = 0 inside the Nyquist 
contour. In this case, the Nyquist criteria requires N = 0 to 
keep the stability for the full system, that is, in the Nyquist 
diagram, the contour of 1 — r s G 0 (£l) should not encircle the 
origin at all. In other words, the contour of r s G 0 (£l ) should 
not encircle the point (1,0) in the (Re[r iS G 0 (Q)],Im[r iS G 0 (Q)]) 
plane. 



